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Abstract 

The internal one-particle density matrix is discussed for Bose-Einstein condensates with finite 
number of particles in a harmonic trap. It is found that the outcome of the diagonalization of 
the single particle density matrix depends on the choice of the internal coordinates: The Pethick- 
Pitaevskii(PP)-type internal density matrix, whose analytical eigenvalues and eigenfunctions are 
evaluated, yields a fragmented condensate, while the Jacobi-type internal density matrix leads to an 
ideal condensate. We give a criterion for the choice of the internal coordinates: In the macroscopic 
limit the internal density matrix should have the same eigenvalues and eigenfunctions as those 
of the corresponding ideal Bose-Einstein condensate in the laboratory frame, this being a very 
physical condition. One choice fulfilling this boundary condition is given by the internal Jacobi 
coordinates, while the internal coordinates with respect to the center of mass of PP do not satisfy 
this condition. Based on our criterion, a general definition of the internal one-particle density 
matrix is presented in a self-bound system, consisting of interacting bosons. 



One of the most amazing phenomena in quantum many-particle systems is Bose-Einstein 
condensation. A characteristic feature of the phenomenon is the macroscopic occupation 



lomogeneous interacting boson 
l| and Yang [2J] who introduced 



of a single-particle state. The criterion of condensation in 
system was given many years ago by Penrose and Onsanger 
the concept of off-diagonal long-range order (ODLRO). The condensate fraction is given by 
diagonalizing the one-particle density matrix which, in general, is defined in the laboratory 
frame. A dominant eigenvalue of the one-particle density matrix implies ODLRO. In other 
words, a system shows Bose-Einstein condensation (BEC) if only one eigenvalue is of the 
order of the particle number N in the system and the others are of the order N° = 1. In 
the case of the existence of several large and comparable eigenvalues, the system is said to 



be a fragmented condensate 



4]. This conventional criterion works well, for example, for 



superconductivity of metals at low temperature and superfluidity of liquid 4 He, including 
recent realization of BEC of ultracold atoms in magnetic traps [5|. 

A controversy concerning the criterion has, however, arisen for the case of condensation 
with attractive interactions among bosonic particles. Wilkin et al. discussed the lowest 
excitation of a condensate with attractive interactions whose system rotates around its 
center of mass but keeps all relative degrees of freedom in its ground state [6|]. It was 
found that the corresponding one-particle density matrix in the laboratory frame has many 
eigenvalues of comparable size, and thus the system should be characterized as a fragmented 
condensate. Pethick and Pitaevskii (PP) claimed that their result is incorrect because only 
the center-of-mass degree of freedom is excited and the internal degrees of freedom remain 
in the lowest states [7|. This is, indeed, physically reasonable. They pointed out that a 
careful treatment is needed in formulating a criterion for Bose-Einstein condensation. They 
proposed that a convenient way to describe correlations in the system is by defining an 
internal one-particle density matrix in which the center-of-mass motion is eliminated, and 
that this definition should give a single eigenvalue of the order of the number of particles. 
Unfortunately, they presented only their idea and did not demonstrate explicit results of 
the diagonalization of their definition of the internal one-particle density matrix. On the 
other hand, Gajda stressed that PP's criterion is not complete in the case of attractive 
systems whose center-of-mass motion has a large uncertainty, namely, whose center-of-mass 
extension significantly exceeds other length scales concerning the internal coordinates 



Recently, a-particle 'condensation' in nuclear systems [2j has attracted much interest. A 



19! 



typical well established example is the 3a condensate in the Hoyle state, i.e. the 0^~ state 
in 12 C, located just above the 3a breakup threshold. It is the finite size counterpart of 
macroscopic a-particle condensation in infinite nuclear matter at low densit y 10f . The 
Hoyle state which was well described by several papers in the past [111 . Il2l . 1 131 ll4j has, 
more recently, been investigated in more detail [9|, [l5|, [ly] and found to have about 1/3 of 
saturation density of normal nuclei and to be, in good approximation, a product state of 
three a-particles, condensed into the lowest mean field 05*-orbit with respect to the motion 
of their center-of-mass coordinates 17|, [l8| . In addition current theoretical work indicates 



that also the Og state of O above the 4a breakup threshold has dilute 4a structure of 



condensate type [19| and possibly states in heavier nuclei as well. 

The character of these a-particle condensates may be seen as a few particle analogue 
to Bose-Einstein condensation of ultracold atoms in magnetic traps [5j. However, some 
qualitatively different features exist between the two systems. Besides the small number 
of particles, there is, for instance, the fact that the a-particles form self-bound systems, 
where the center-of-mass motion can not be controlled by the external field, because the 
total wave function of the system should, in principle, be described in terms of only the 
internal coordinates, eliminating the center-of-mass degree of freedom. Thus, the issue of 
the internal one-particle density matrix with a definite finite number of particles is very 
relevant in the study of a-particle condensation in finite nuclear systems. 

The purpose of this paper is to give a suitable definition of the internal one-particle den- 
sity matrix of a self-bound system with a finite particle number. The criterion by PP seems 
to work well for this issue at first sight. However, we will demonstrate that their criterion 
is not adequate, leading to a fragmented condensate, contrary to their initial objective. For 
illustration, we treat a simple case, i.e. the internal state of a Bose-Einstein condensate with 
finite particle number in a harmonic trap. Let us first consider the one-particle density 
matrix in the laboratory frame for an ideal Bose-Einstein condensate with N spinless parti- 
cles. The result is trivial but instructive for studying the nature of the internal one-particle 
density matrix, as will be discussed later. 

The iV-particle Hamiltonian in laboratory frame is presented as 

H^Pl + E^rl (1) 

1=1 1=1 

The ground-state wave function of this system is expressed as a product of identical Gaussian 



A /2zA 3 / 4 



single-particle wave functions, i.e. 

*({r}£i) = n(-)'^xp(-,r 4 2 ), (2) 

where v = muj/2h, and {r}^ denotes the set of the coordinates r, (i = 1 • • • , AT). The one 
particle density matrix in the laboratory frame is defined as 

N 



p£kr,r') = /n^«**(r,{r«}£ 3 )*(r / > {r i }^ 2 ), (3) 

J i=2 

'2v\ 



4 = 2 

3/2 



7T 



■ 



exp 



-z/(r 2 + r' 2 ' 



(4) 



It is noted that the density matrix is independent of the number of particles N and is 
separable with respect to r and r'. The separability originates from the fact that the 
Hamiltonian is separable for r j and p t in Eq. ([TJ . 

The nature of the single particle orbits and their occupation probabilities in the relevant 
system can be obtained by solving the eigenvalue problem of the density matrix (j4]). This 
is easily done, and we find that the density matrix has only one non zero eigenvalue with 
one eigenf unction, namely, the zero-node S"-wave Gaussian y?(r*) = (2z//7r) 3//4 exp(— ur 2 ) (or 
OS* harmonic oscillator wave function R<m{r,v)) with 100 % occupancy. This means that 
all particles are condensed in that single orbit, i.e. an ideal Bose-Einstein condensation is 
realized in the laboratory frame. To say it again, this feature is independent of the number 
of particles N. 

Next we consider the internal one-particle density matrix for the iV-particle Bose-Einstein 
condensation in a harmonic trap described by the wave function Eq. (J2J) with the total Hamil- 
tonian Eq. (CQ) . Internal means that the density is free from the center-of-mass coordinate 
of the system. In the present paper, two kinds of internal coordinate sets are introduced, 
1) coordinates with respect to the center of mass and 2) Jacobi coordinates. The former set 
was first considered by Pethick and Pitaevskii [7] to define the internal one-particle density 
matrix. We call it Pethick-Pitaevskii-type (PP-type) internal one-particle density matrix in 
the present paper. For the latter set, we call it Jacobi- type density matrix. 

PP-type internal one-particle density matrix 

In order to define an internal one-particle density matrix, Pethick and Pitaevskii adopted 
internal coordinates defined with respect to the center of mass [TJ. The center-of-mass 



coordinate R and the coordinate q i of particle i relative to the center of mass are given, 
respectively, by 



1 



N 



R =jjY r ^ Qi = ri-R, 

i=l 



(5) 



where the definition of the center-of-mass coordinate R implies that only N — 1 of the q i 
are independent. Here, we define the conjugate momenta 7Tj and P for the coordinates q t 
and R, respectively. Then, the total Hamiltonian in Eq. (TjQ) is rewritten as 



H = H mt + H t 



(6) 



2m 
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N 

N-l 



N-l 



N-l 
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i=l iv i<i'=l 



H r 



N-l 

<'<''-" Y Qi 2 + Y Qi-Qv 

i=l i<i'=l 

1 P 2 + -Nmu 2 R 2 , 



2Nm ' 2 

where H mt and H cm denote the internal and center-of-mass Hamiltonians, respectively. 
The total wave function in Eq. (J2J) is expressed as 

1 



(7) 
(8) 



*({ri}£i) 



AT3/2 



$in t ({gjf=l) 



X ®int({qi}?=l l )$ C rn(R) 

' N x (2u) N - 1 \ 3/4 



(9) 



7T 



JV-1 



x exp 

®r. m (R) 



N-l 



- Y i 5 H' + l ) V( ii ■ Qi 1 



i,i'=l 



(^feM-N»R% 



(10) 
(11) 



where $i nt and $ cm denote the internal and center-of-mass wave functions, respectively. 
According to Pethick and Pitaevskii |7| , the internal one-particle density matrix is given 

as y, 



Pinippte' q') = J dq 2 --- dq N _ 1 



N \ 3 / 2 /2u\ 3/2 
N-l) V W 



x exp 



3iV-2 

~2(N - 1) 



V\ 



HQ + Q ) + 



N -2 
N-l 



vqq 



(12) 



(13) 



It is noted that this density matrix depends on the number of particles N and contains the 
cross term q ■ q'. The origin of the cross term comes from the nonseparability of 7Tj and q L 
in the internal Hamiltonian [see Eq. (|7j)]. One shall remark that the result of Eq. (1131) differs 
substantially from Eq. (9) in the paper by Zinner and Jensen 21]. 

Let us discuss the nature of the internal one-particle density matrix pl n l(q,q'). First 
we study the single-particle orbits and their eigenvalues obtained by solving the eigenvalue 
problem, 

J(^(q,4)<p(qW = Mq)- ( 14 ) 

We find that this equation can be solved analytically. The single-particle orbits tp are 
expressed by the harmonic oscillator wave functions RnLM^Q, Pn) with the orbital angular 
momentum L, magnetic quantum number M, quanta Q = 2n + L (n = 0, 1, • • ■), and size 



parameter j3^ = J2N/(N — l)v. The eigenvalues or occupation probabilities A are given 
as [22 1 



x (lm) = (4Nf/\N-2)^ L 

n,N [3N-2 + 2^2N{N - \)fn+L+z/2 ' 



and satisfy the following completeness relation, 

oo L oo 

E E EC f) = i. (i6) 

L=0 M=-L n=0 

The occupation probability with respect to the partial wave with quantum number L is 
defined as 

L oo 

Ai?= E E« f) - (17) 



M=-Ln=0 



In the macroscopic limit 



ASv L L = (2L + 1)(2 - v^)(3 - 2v/2) L . (18) 



We remark that the summed eigenvalues A^^ still depend on the angular momentum L 
(see Fig. HJ. 

Jacobi-type internal one-particle density matrix 




FIG. 1: Spectra of the occupation probabilities A^ for the PP-type internal density matrix [see 
Eq. (HID]. 



For the iV-particle system, we define the N — 1 internal Jacobi coordinates {£j}j = i and 
the center-of-mass coordinate R as follows: 



1 



N-i 



N 



rN-i+i 



N-i 



5>fc> r =mJ2 



(19) 



fe=l ^ i=l 

where £± denotes the relative coordinate between the N-th particle and the remaining (N—l) 
particles, and other Jacobi coordinates are self-evident. Then, the iV-particle Hamiltonian in 
Eq. flTJ) can be separated into the internal and center-of-mass Hamiltonians, H = H mt + H cm , 
where 

"^ 7T? , »=} W*g N-i 

t*v* ■ h 2 >*- N -i + r 



H; 



int. 



-m, 



(20) 



7Tj denotes the conjugate momenta with respect to coordinate £ i5 and H cm is given in Ec^ 



The Jacobi-type one-particle density matrix with respect to ^ l and £[ is given as |20| 



&&0 =/<*&• ••#*-! 



N-I2v\ 3/2 

1 exp 



iV 7T 



■ v -V + «*: 



iV 



(21) 
(22) 



where $mt,j represents the internal, fully symmetric [13, [l_8|, |20|], wave function in Jacobi 
coordinates corresponding to Eq. (J2J). This choice of the coordinate ^ 1 for the internal 
density matrix is natural, because the single particle orbit should be defined with respect 
to the relative coordinate between one particle and the other remaining N — 1 particles in 



the Jacobi coordinate system. The separability with respect to $, and £' in Eq (|22|) stems 
from the fact that Jacobi coordinates form an orthogonal system (it is noted that the cross 
term q ■ q' does appear in the case of PP [see (TT31) ] because the internal coordinates used 
by PP are not orthogonal). The eigenvalue equation of the one-particle density matrix 
can be solved analytically. We find that the density matrix has only one eigenfunction 
<Pm,J,N = Rooo(£,, (N — l)u/N) with non-zero eigenvalue, namely, the OS harmonic oscillator 
wave function with 100 % occupancy, that is, in the same notation as above 

Aiv L) = E ^ LM) = ho- (23) 

M=-L 

This means that all particles are condensed in the single OS" particle state independent of N, 
although the size parameter in the eigenstate <fwt,j,N depends on N and is slightly different 
from that in the eigenfunction Rooo{r, v ) m the laboratory frame, as discussed above. 

We, therefore, see that the eigenvalues and eigenfunctions of the internal density matrix 
depend on the choice of the internal coordinates. This is a somewhat surprising result. In 
order to overcome the difficulty, we give a criterion for the choice of the internal coordinates: 
In the macroscopic limit (N —>■ oo) the internal density matrix should have the same eigen- 
values and eigenfunctions as those of the corresponding ideal Bose-Einstein condensate in 
the laboratory frame. This is a very physical boundary condition. Obviously, the PP-type 
one-body density matrix does not satisfy the condition, while the density matrix of the 
Jacobi-type does. These results mean that one should take internal coordinates which do 
not produce any correlation in the internal one-particle density matrix in the macroscopic 
limit. Otherwise, unphysical situations occur like the density matrix of the PP-type which 
clearly is inadequate. 

The present considerations can be applied to a general case of the internal one-particle 
density matrix for interacting bosons in a self-bound system. The internal Hamiltonian of 
the system or translationally invariant shell-model Hamiltonian with the Jacobi coordinates 
(TT9|) is presented as H int = H — H crn , 

N-l 2 N-l .. ,.2/2 N 

tf- = E j- + E ^T~ + E «(r, - r,), (24) 

where the center-of-mass Hamiltonian H cm of Eq. (jHJ) is subtracted, and v presents the 
residual two-body interaction between bosons. Defining ^intd^}^! 1 ) as the eigenfunction 
of H int , the internal one-particle density matrix of the system is presented as 



-,(!> 



pffi&O 



JV-1 



n<% 



i=2 
I JV-1 



iV-l^ 



x*k(M*J«)<Me',te}£ 



(25) 



which is similar to Eq. ( 1221) . This definition satisfies the physical boundary condition in the 
limit of N —>■ oo in the case of the absence of the interactions among bosons. Full symmetry 
of the internal one-particle density matrix (1231) can be verified 17|, [18|, [23]. 

In conclusion, we discussed the internal one-particle density matrix for Bose conden- 
sates in a harmonic trap and for self-bound Bose-systems. It was found that the PP-type 
one-particle density matrix is physically inadequate to study the internal degree of Bose 
condensation, while the density matrix of the Jacobi-type is fully appropriate. The use of 
the latter is, thus, of great importance for the exploration of condensates in self-bound sys- 
tems such as a-particle condensates in nuclei or small droplets of superfluid 4 He, where only 
internal degrees of freedom are relevant. Indeed, the definition (123|) has successfully been 
applied to study the degree of Bose condensation for self-bound a- particles in nuclei, un- 
ambiguously demonstrating that self-conjugate nuclei such as 12 C, 16 0, ■ ■ • show long-lived 



excited states close to the na disintegration threshold where the na 



approximation, a product state of OS-orbitals, that is a condensate |l7l . 



partic 



18 



cs 



brm, in good 
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